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We study the character of the first excited l/2 + state of 9 Be, which is observed as a low-lying 
sharp peak in the cross section of 9 Be(7,n)2a just above the s Be + n threshold. Using the a + a + n 
three-body model, we describe the ground and excited unbound states of 9 Be above the a + a + n 
threshold. Applying the complex scaling method to the three-body model, we find no l/2 + resonant 
solutions with the scaling angle of 6 < 15 degrees, while the low-lying peak in the photodisintegra¬ 
tion cross section is reproduced in the present calculation. It is found that the low-lying peak is 
dominantly explained by the s Be + n component. Furthermore, using the analytical continuation 
of the coupling constant of the three-body interaction for the a + a + n system, we discuss the 
virtual-state character of the l/2 + state. 

PACS numbers: 21.45.-v, 21.60.Gx, 26.30.-k, 27.20.+n 


Introduction. The neutron capture (n, 7) and its in¬ 
verse (7, n) reactions are interesting topics in nuclear 
physics since these reactions play important roles in nu¬ 
clear astrophysics. In neutron capture reactions, the 
a(an, 7) D Be reaction is one of the most interesting ones. 
It has been suggested that the reaction rate of this re¬ 
action is crucial to understand the productions of heavy 
elements in supernova explosions [l], 0 • 

To determine the reaction rate of the a(an, 7) 9 Be, the 
9 Be(7,n) 8 Be reaction has been performed by several ex¬ 
periments 0-0. In the a(an, 7) 9 Be reaction, a sequen¬ 
tial process, a(a, 7) 8 Be(n, 7) 9 Be, has been considered as 
a dominant one. However, owing to the short life-time 
of the 8 Be ground state (~ 10 _16 s), a direct measure¬ 
ment of the 8 Be(n, 7) 9 Be reaction is impossible. For an 
alternative way, the cross section of its inverse reaction, 
9 Be(7, n) 8 Be, has been measured to deduce the cross sec¬ 
tion of 8 Be(n, 7) 9 Be. 

In the experiments, there is an inconsistency among 
the observed cross sections of 8 Be(n, 7) 9 Be in particular 
the peak just above the 8 Be + n threshold 0-0. The 
main difference among the cross sections comes from the 
contribution from the first excited l/2 + state of 9 Be via 
the E 1 transition. This low-lying l/2 + state can have an 
impact on the reaction rate of 8 Be(n, 7) 9 Be in stellar en¬ 
vironments and supernova explosions [l], 0 ■ For theoret¬ 
ical side, it is interesting problem to answer how the low- 
lying l/2 + state of 9 Be contributes to the 8 Be(n,7) 9 Be 
reaction. 

Theoretically, the structure and the reaction mecha¬ 
nism of 9 Be have been studied by using the a + a + n 
three-body models 00. In previous works, the struc¬ 
ture of the first excited 1/2+ state has been discussed. 
Efros et al. 0 discussed the structure of l/2 + state from 
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the scattering length for the 8 Be + n scattering. They 
show that the scattering length is obtained as a large 
negative value, which characterizes the first excited l/2 + 
state as a virtual state. Also, Arai et al. 0 show that 
the virtual-state character of the l/2 + state from the R- 
matrix analysis. On the other hand, in Refs. [III. IT0 . 
the first excited l/2 + state is discussed as a three-body 
resonance of a + a + n. In these calculations, they 
show the importance of a strong mixture of the 5 He + a 
configuration at the internal region of the l/2 + state. 
There exists the contradictions among the theories and 
the structure of the first l/2 + is still unclear. Addition¬ 
ally, in Ref. 00, the importance of the direct three- 
body capture for 9 Be below the 8 Be + n threshold is sug¬ 
gested recently. It is important to investigate the relation 
between the structure of the first excited 1/2+ state and 
the mechanism of the capture reaction. Comprehensive 
understanding for the structure of l/2 + of 9 Be is now 
required. 

The purpose of this work is to understand the structure 
of the low-lying l/2 + state of 9 Be and its contribution 
to the 9 Be(7, n) 8 Be reaction. We describe 9 Be by using 
the a + a + n three-body model. To treat the pho¬ 
todisintegration into unbound excited states of 9 Be, we 
apply the complex scaling method (CSM) 0-0 to the 
a + a + n three-body model. The CSM is a powerful 
tool to investigate many-body resonances and has been 
extensively used to discuss their structures 0, 0, 0 . 
It has been shown that the CSM enables us to describe 
the resonances from two- to five-body system success¬ 
fully 0 - 0 | . Furthermore, the CSM has been developed 
to describe the breakups of unstable nuclei. The CSM has 
been applied to the Coulomb breakup reactions of two- 
neutron halo nuclei into core + n + n scatterin g st ates 
and well reproduces the observed cross sections 0-0]. 
From the above, the CSM is promising to discuss the 
excited states of 9 Be and the photodisintegration into 



2 


unbound excited states of * 9 Be on the same footing. 

Generally, the resonances are obtained as the pole of 
5-matrix with complex energy eigenvalues in the fourth 
quadrant of the complex energy plane. On the other 
hand, the virtual state is obtained on the negative en¬ 
ergy axis of the second Riemann sheet. In the CSM, 
we can obtain the solution of the resonances but not 
of the virtual state. Such a pole of the virtual state 
has been discussed by using an analytical continuation 
method within the two-body system [101 ]. The properties 
of the virtual state of a three-body system such as 9 Be 
is still unclear. In this work, we investigate not only the 
eigenvalue distribution but also the photodisintegration 
cross section in order to clarify the structure of the l/2 + 
state whether the resonance or the virtual state. 

Formalism. We solve the Schrodingier equation for the 
a + a + n system using the complex-scaled orthogonality 
condition model [13 • The complex-scaled Schrdingier 
equation is given as 


charge radius of the a particle. In the present calculation, 
A is taken as 10 6 MeV. 

To discuss the photodisintegration of 9 Be, it is im¬ 
portant to reproduce the breakup threshold into the 
a + a + n. In the present calculation, we introduces 
the a + a + n three-body potential V 3 to reproduce the 
binding energy of the 9 Be ground state, E gs , measured 
from the a + a + n threshold. The explicit form of V 3 
is given as 

V 3 = u 3 exp(-/qo 2 ), (6) 

where p is the hyper-radius of the a + a + n system. 

We solve the eigenvalue problem given in Eq. 0 with 
Gaussian expansion method [3lj| - and obtain the energy 
eigenvalues and eigenstates (their biorthogonal states) 
as {E 8 } and ({'Pj($)}), respectively. Using 

them, we define the complex-scaled Green’s function 
Q 6 {E-i,i') as 


H e ^j(9) = E e v Wj{9), (1) 

where J is the total spin of the a + a + n system and v 
is the state index. The complex-scaled Hamiltonian and 
wave function are given as 

H 8 = U(9)HU~ 1 (9) and \E^(6>) = U(0)* v j, (2) 

respectively. The complex scaling operator U(9) trans¬ 

forms the relative coordinate £ as 

U(9) (3) 

where 9 is the scaling angle being a positive real number. 

The Hamiltonian for the relative motion of the 
a + a + n three-body system for 9 Be is given as 

3 2 

H = T c . m . + ^ V an {&) + V aa + U PF + V 3 , (4) 

i=l i=1 

where t-i and T c m . are kinetic operators for each particle 
and the center-of-mass of the system, respectively. The 
interactions between the neutron and the i-th a particle 
is given as V an (£i), where is the relative coordinate 
between them. We here employ the KKNN potential (28j 
for V an ■ For the a-a interaction V aa we employ a folding 
potential of the effective NN interaction (2!| and the 
Coulomb interaction: 


S e (£;€,£') 


1 

E-H 9 




^ v {9)^ v {9) 

E-El 


(7) 


In the derivation of the right-hand side of Eq. 0 , we 
use the extended completeness relation, whose detailed 
explanation is given in Ref. (32}. It is noticed that we 
take into account outgoing boundary conditions for all 
open channels of a three-body system in the form of com¬ 
plex energy eigenvalues E 8 . The complex-scaled Green’s 
function in Eq. 0 enables us to describe the scattering 
observables for many-body systems, such as the photo¬ 
disintegration cross section. 

We calculate the cross section of the photodisintegra¬ 
tion of 9 Be(3/2 _ ) + 7 — > a + a + n in terms of the 
multipole response. The cross section is expressed as the 
following form; 




(2tt) 2 (A + 1) (E^- 1 dB(E\,E 7 ) 

A[(2A + 1)!!] 2 \Hc) dE 1 ’ U 


where E 1 is the incident photon energy and B(E\) is 
electric transition strength with the rank A. We here 
calculate the photodisintegration cross section from the 
ground 3/2“ state to l/2 + states in 9 Be and consider 
only the El transition here. Using the CSM and the 
complex-scaled Green’s function in Eq. 0, the El tran¬ 
sition strength is given as 


V aa (r) = exp (-ar 2 ) + —erf (fir ), (5) 

where vq = —106.09 MeV, a = 0.2009 fm -2 , and f3 = 
0.5972 fm -1 . The pseudo potential Vpf = A|$pf)( < 1 ) pf| 
is the projection operator to remove the Pauli forbidden 
states from the relative motions of a-a and a-n [3Cj. The 
Pauli forbidden state is defined as the harmonic oscilla¬ 
tor wave functions by assuming the (0s) 4 configuration 
whose oscillator length is fixed to reproduce the observed 


X E — E 8 (^ 4 / 2+ (0)||O e (El)||vl> gs ), 

(9) 

where J gs and Tgs represent the total spin and the wave 
function of the ground state, respectively. The energy E 
is related to E^ as E = E^ — E gs . From the Eqs. 0 
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TABLE I. Ground-state properties in comparison with exper¬ 
iments. The calculated binding energies (E gs , unit in MeV), 
charge radii (R c h, unit in fm), and matter ones (R m , unit in 
fm) with and without the three-body potential V3 are listed. 



Eg a 

Rch 

Rm 

V3 

d 

with V3 

1.57 

2.53 

2.42 

1.10 

0.02 

without V3 

2.14 

2.50 

2.39 

- 

- 

Exp. 

1.57361 

2.519±0.012 1 

2.38T0.011 




a Reference 
b Reference 
c Reference 


m 

34] 

m 



and ©, we finally obtain the photodisintegration cross 
section as 




16tt 3 ^ dBiEl^E^) 
9 He 7 dE 1 


( 10 ) 


Results. We first show the calculated ground-state 
properties of 9 Be. In TABLE m the calculated binding 
energy and charge and matter radii are listed. Without 
the three-body potential in Eq. ©, the binding energy 
of 9 Be ground state is overbound and the charge radius 
is slightly small compared to experiments. To reproduce 
these quantities, we need the repulsive three-body poten¬ 
tial whose parameters are given as V 3 = 1.10 MeV and 
/x = 0.02 fnG 2 . As results, we reproduce the binding 
energy and charge radius of 9 Be ground state simulta¬ 
neously, while the matter radius is slightly larger than 
the observed one. We also confirm that no resonance of 
the 1/2+ state is found with the three-body potential in 
TABLE Q] for the 9 = 15 degrees case. 

Next we discuss the photodisintegration cross section 
to l/2 + states. In the present calculation, we fix the 
ground-state wave function obtained with the three-body 
potential in TABLE [U In FIG. [l] we show the calculated 
cross sections using Eq. (11011 in comparison with the two 
sets of the observed data which commonly have peaks 
just above the 8 Be + n threshold. The red (dashed) 
and blue (dotted) lines show the cross sections with and 
without the three-body potential for excited l/2 + states, 
respectively, whose parameters are same as those in TA¬ 
BLE HI In both results, the calculated cross sections un¬ 
derestimate the low-lying peak above the 8 Be + n thresh¬ 
old. 

To discuss the observed sharp peak just above the 
8 Be + n threshold in the photodisintegration cross sec¬ 
tion, we change the strength, V3, for the l/2 + state to fit 
the observed data but its range is fixed as the same as 
used in the ground state. We here take the strength as 
V 3 = —1.02 MeV for the l/2 + state and obtain the cross 
section as shown as the black (solid) line in FIG. [T] Our 
result reproduces the observed peak by using the attrac¬ 
tive three-body potential. We confirm that the calculated 
cross section rapidly increases just above the 8 Be + n 
threshold and there is negligibly small strength below 
this threshold. We also find that the calculated cross 


FIG. 1. Calculated photodisintegration cross sections in com¬ 
parison with the experimental data. The red (dashed) and 
blue (dotted) lines are results with and without the three- 
body potential, where V 3 = 1.10 MeV and /x = 0.02 fm. The 
black (solid) line represents the cross section calculated by us¬ 
ing an attractive three-body potential with V 3 = —1.02 MeV. 
The experimental data below E 1 = 2.2 MeV are taken from 
Refs. [f| and @]. The arrow indicates the threshold energy of 
the 8 Be(0 + ) + n channel. 



FIG. 2. Distribution of energy eigenvalues of l/2 + states 
measured from the a + a + n threshold with scaling angle 
9 = 15 degrees. We here employ the three-body potential with 
V3 = —1.02 MeV and fi = 0.02 fm~ 2 . The red (solid), blue 
(dashed), and green (dotted) lines represent the branch cuts 
for a + a + n, 8 Be(0 + ) + n, and s He(3/2 _ ) + a continua, 
respectively. 


sections show the strong dependence on the strengths of 
the three-body potentials as shown in FIG.[H This result 
is interesting and suggests the existence of the three-body 
unbound state of 9 Be(l/2 + ), such as a resonance or vir¬ 
tual state. In relation to the cross section, we discuss the 
character of the l/2 + state. 

We investigate the origin of the low-lying peak above 
the 8 Be + n threshold more in detial. For this purpose, 
we show the distribution of the energy eigenvalues of the 
l/2 + states by using the CSM. In the CSM, continuum 
states are obtained along the branch cuts which start 
from the threshold energies and are rotated down by 29. 
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FIG. 3. Decomposed photodisintegration cross sections. The 
red (solid) and blue (dashed) lines are contributions of the 
a + a + n and s Be + n continuum states. The black thin 
line is same as that in FIG. [T] 


A resonance is obtained as a solution with complex en¬ 
ergy of E e = E r — iT/2 isolated from the continuum ones 
. On the other hand, the virtual states and broad reso¬ 
nances, which are located on the second Riemann sheet 
covered by the rotated first Riemann sheet, cannot be 
obtained as the isolated pole in CSM. The contributions 
from these states to the cross section are scattered into 
the continuum states rotated on the 29 lines. In FIG. [2j 
we show the distribution of the energy eigenvalues for 
the l/2 + states calculated with V 3 = —1.02 MeV, which 
reproduces the observed peak as shown in FIG. [D In 
the present calculation, we find no resonances in the en¬ 
ergy eigenvalue distribution. All energy eigenvalues are 
located on the 20 -lines, corresponding to the branch cuts 
for the a + a + n, 8 Be(0 + ) + n, and 5 He(3/2 _ ) + a 
continuum states. 

We investigate the contributions of two- and three- 
body continuum states to the cross section to under¬ 
stand the mechanism of the photodisintegration. Using 
the CSM, we decompose the cross section into 8 Be + n 
and a + a +n components as shown in FIG. [3j to clar¬ 
ify the dominant contribution in the cross section |24| . 
The 5 He + a contribution is found to be negligible in 
the low-lying region, and we do not show it in FIG. [3] 
From the results, we see that the 8 Be + n component 
is almost identical to the total cross section. This fact 
indicates that the 8 Be + n decay is dominant in the pho- 
todintegration. This decay mode should be related to the 
structure of the 1/2+ state of 9 Be. 

To investigate the structure of the l/2 + state of 9 Be, 
we calculate the energy eigenvalues of the a + a + n sys¬ 
tem by changing the strength of the three-body potential, 
V 3 , which is shown in FIG.|4] In the present calculation, 
when the strength of the three-body potential V 3 = —1.3 
MeV, the resonance pole suddenly appears just below 
the 8 Be(0 + ) + n threshold. This resonance pole with a 
narrow decay width moves smoothly to the bound state 
region as the three-body potential becomes more attrac¬ 
tive, and we finally obtain the 9 Be bound state with the 



FIG. 4. Pole trajectory of the 9 Be l/2 + state in complex en¬ 
ergy plane by changing the three-body potential. The closed 
circles represent the poles obtained as isolated three-body res¬ 
onances in CSM. The open circles and the crosses are spec¬ 
ulated pole positions for the virtual states and broad reso¬ 
nances, respectively. 


region of V 3 < —1.8 MeV. On the other hand, we consider 
the pole trajectory in the opposite case of the three-body 
potential with V 3 > —1.3 MeV. If the resonance exists, 
the pole with a narrow decay width should appear above 
the 8 Be + n threshold as the analytical continuation from 
the resonance pole as shown with the crosses in FIG. [4] 
However, we found that no resonances appear above the 
8 Be(0 + ) + n threshold for V 3 > —1.3 MeV of the three- 
body potential. These facts in the pole trajectory show 
the possibility of the virtual state of the l/2 + state con¬ 
sisting of 8 Be(0 + ) + n when we take V 3 = —1.02 MeV, 
which reproduces the experimental cross section. The 
existence of the virtual state is consistent with the dom¬ 
inant decay into 8 Be + n in the photodisintegration of 
9 Be. 

The virtual state is often discussed in unstable nuclei, 
such as 10 Li and 11 Li [25], [2(1 [36]. For 10 Li, Masui et 
al. j3(| showed the explicit pole trajectory of the virtual 
state of 10 Li as 9 Li + n picture, which is similar to the 
present 8 Be(0 + ) + n situation. In FIG[H we schemati¬ 
cally plot the pole trajectory in the region of V 3 > —1.3 
MeV, considering the analogy with the 10 Li case as shown 
in FIG. 3(b) of [36j]. This trajectory is located on the sec¬ 
ond Riemann sheet of the 8 Be + n system which cannot 
be obtained in the CSM. 

We consider the different case when there exists the 
l/2 + resonance to see how the resonance contributes 
to the photodisintegaration cross section. By using the 
three-body potential with V 3 = —1.3 MeV, the 1/2+ res¬ 
onance is obtained with the energy and decay width be¬ 
ing 0.091 and 0.002 MeV, respectively, just below the 
8 Be + n threshold at 0.0918 MeV measured from the 
a + a + n threshold. In FIG. [U we show the cross sec¬ 
tion in this case, which shows quite a sharp peak at the 
resonance energy. From the decomposition of the cross 
section, we find that the peak in the cross section is dom¬ 
inated by the resonance while the 8 Be + n continuum 
states have a sizable contribution to the cross section. 
This trend is much different from that in FIG. [3] In the 
result of FIG. [3j the peak in the cross section is con- 
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FIG. 5. Decomposed photodisintegration cross sections by us¬ 
ing the three-body potential with V 3 = —1.3. The red (solid), 
blue (dashed) lines are contributions of the a + a + n and 
8 Be + n continuum states. The green (dotted) line is that of 
the resonance pole. The black thin line represents the sum 
of them. The arrow indicates the threshold energy of the 
8 Be(0 + ) + n channel. 


structed by several low-lying eigenstates located on the 
8 Be + n continuum 28 line in FIG. [2] 

Summary. We investigate the character of the 1/2+ 
state of 9 Be using the photodisintegration reaction with 
the a + a + n three-body model and the CSM. The cal¬ 
culated photodisintegration cross sections into the l/2 + 
states are shown to have a strong dependence on the 
strength of the three-body potential for the 1/2+ state. 
The experimental cross section shows the sharp peak 
just above the 8 Be + n threshold, which is nicely re¬ 
produced with the attractive three-body potential. We 
cannot find any resonance poles for the l/2 + states in 
explaining the peak in the cross section. From the de¬ 
composition of the calculated cross section, it is shown 
that the 8 Be + n continuum states dominates the cross 
section to the l/2 + states. These results indicates the 
possibility of the virtual-state nature of the first excited 
1/2+ state. In addition, the pole trajectory suggests that 
the pole of the l/2 + state is located on the second Rie- 
mann sheet of 8 Be + n instead of the broad resonances 
of a + a + n. 
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